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ABSTRACT 



The vertical distribution of ozone was observed by the Air Force 
Cambridge Research Laboratories North American ozonesonde network from 

29 April 1963 to 10 May 1963 on a daily basis. This observed data, 
supplemented by twenty-four hour trajectory calculations, was used 

to prepare distributions of ozone mixing ratio on the 200, 100, 50, and 

30 mb surfaces. The ozone distributions were evaluated for the ozone 
transport by mean and eddy motions in the lower stratosphere. The 
results of these calculations show the magnitude of the horizontal 
and vertical transport of ozone into and out of the 3^/50, 50/100, 
and 100/200 mb volume over North America. The mean vertical trans- 
port, supported by the vertical eddy transport, is combined to move 
the ozone primarily downward through the stratosphere. At the tropo- 
pause the vertical eddies transport the ozone to the troposphere for 
eventual destruction. 
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I, INTRODUCTION 



Ozone is an important but not well-documented atmospheric gas. 

Its importance to meteorologists is twofold. First, ozone can be 
used as a quasi -conserva tive atmospheric trace substance below 30 
kilometers and, second, as an absorber of ultraviolet light, ozone 
affects the heat balance through the diabatic heating term of the 
first law of thermodynamics. Although the exchange of energy by 
ozone is worthy of investigation, ozone as an atmospheric tracer is 
the subject of this study. Before proceeding, a review of the 
significant works of earlier investigators is desirable. 

The first major investigation of the measurement and distribution 
of atmospheric ozone was carried out under the direction of G. M. B. 
Dobson during a five-year period from 1924 to 1929. Dobson and his 
associates reported the findings of their study in a series of four 
articles in the Proceedings of the Royal Society of London (1926, 
1927, 1929, 1930). Based on the work of previous investigators, 
Dobson developed a spectroscopic instrument for use on the earth's 
surface that measures the total amount of atmospheric ozone in a 
vertical column. The results of daily measurements of total ozone 
at Oxford during 1924 and 1925 were compared with surface pressure. 

It was found that surface pressure had a high negative correlation, 
on the order of -.5, with total ozone measurements. On an annual 
basis, a maximum of total ozone was observed in the spring with a 
corresponding minimum in the fall. Comparisons of total ozone were 
made with atmospheric pressure at nine, twelve, and fourteen 
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kilometers. A high negative correlation, on the order of -.8, was 
observed . 

In 1926 instruments were sent to various locations in north- 
western Europe. The observations produced similar conclusions as 
the Oxford findings of 1925. Using the results of 1926-1927 a com- 
plete description of the ozone conditions relative to cyclones and 
anticyclones was presented. In general, high total ozone was measured 
in the area of northerly wind flow and low total ozone in the sectors 
associated with southerly wind flow. Highest total ozone values were 
found to the southwest of the surface low pressure centers. 

During 1928 and 1929 the spectrographs were redistributed so that 
observations were taken at Arosa (Switzerland), Table Mountain 
(California), Helwan (Egypt), Kodaikanal (India), Christchurch (New 
Zealand), Oxford (England), and Montezuma (Chile). The analysis of 
the measurements showed that the midlatitude stations had similar 
observations as discussed for northwestern Europe. The tropical 
stations observed no significant pressure or total ozone variations. 

Several years later, using the observations obtained during the 
period 1928-1932 in northwestern Europe, Meetham (1937) found a high 
positive correlation between total ozone and temperature at a height 
of 12, 15, and 18 kilometers. The magnitude of the correlation 
coefficients was on the order of +.5. He also found a high negative 
correlation between total ozone and pressure at 9, 12, 15, and 18 
kilometers altitude and between total ozone and the height of the 
tropopause. Although values vary somewhat, the pressure correlation 
coe f f i c i ent s have a magnitude of -.4, while the corresponding tropo- 
pause figures are -.5. 
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Godson (1960) has summarized Northern Hemisphere total ozone obser- 
vations over many years by month and latitude. His results show a 
strong middle and north latitude spring maximum during the month of 
March and a middle and north latitude minimum during the month of 
October. Relatively constant values over the seasons are shown in 
the tropics with increasing values northward. The maximum value is 
reported above 70N during the month of March. The steepest latitudinal 
gradient occurs across middle latitudes during late winter and early 
spring. The flattest gradient occurs during late summer and early 
autumn. This elaborates upon and confirms earlier and limited obser- 
vations. Godson*s summary is generally accepted as a reasonable 
description of the mean distribution of total ozone. 

The vertical distribution of ozone has been difficult to obtain. 
This is attributed to the lack of a reliably accurate, inexpensive, 
and expendable instrument. A summation of the various methods of 
measuring the vertical distribution of ozone is presented by Craig 
(1965) and therefore is not presented here. The most widely used 
technique has been the Umkehr method devised by F. P. W. Gotz in 
1931 (Craig, 1965). Using photometery, this method breaks the 
atmosphere into five or nine layers. For each layer a single ozone 
value is computed. Of more recent use are the Regener chemiluminescent 
ozonesonde and the Brewer-Mast ozonesonde. Both of these instruments 
are capable of being used on radiosondes and are able to give a much 
finer vertical ozone structure than the Umkehr method. 

Investigators (see, for example, Gotz, 1951; Ramanathan and 
Kulkarni, 1960; Kulkarni, 1962; Mateer and Godson, 1960; Craig, 

1960, 1965; Hering and Borden, 1965; or Bojkov, 1965, 1969) agree on 
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the general shape of the vertical ozone profile. Low, relatively 
constant, mixing ratios of ozone are found in the troposphere with 
a very sharp increase at the tropopause. The maximum value of ozone 
mixing ratio is most often found in the 15 to 30 kilometer layer with 
the altitude of the maximum decreasing poleward. Above this maximum 
the ozone mixing ratio slowly decreases or remains relatively constant 
with altitude. Maximum hemispheric values are observed near the 30- 
kilometer altitude of the equatorial latitudes. In the stratosphere 
isolines of ozone distribution normally slope poleward and downward 
from the tropical regions. Large daily fluctuations in ozone content 
are observed in the region between the tropopause and the primary ozone 
maximum. The fluctuations are weakest in tropical regions and strongest 
in middle and north latitudes. 

It is worthwhile at this point to consider the expected distribution 
of ozone assuming photochemical equilibrium. Gotz (1951), Craig (1950, 
1965), Prabhakara (1963), and Brewer and Wilson (1968) have discussed 
the origin and details of photochemical theory with respect to atmos- 
pheric ozone. Therefore, photochemical theory will not be reviewed 
here. From the work of the previously noted authors, two salient 
points a re noted : 

1. Photochemical equilibrium theory correctly predicts the 
observed amount of ozone in the atmosphere above approximately 30 
kilometers. Predicted by this theory is the equatorial maximum region 
and decreasing values poleward and with altitude. Below this region, 
phot ochemica 1 theory does not agree with observed distributions. 

2. Production and recombination can be considered small below 
the 30 kilometer maximum observed values. Therefore in this region 
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ozone must be quasi-conservative and can be used as an atmospheric 
tracer. 

Since the observed distribution and photochemical theory disagree, 
the discrepancy must be due to transport processes. Transport mechanisms 
must move the ozone from high equatorial source regions to the lower 
stratosphere in middle and high latitudes. From the lower stratosphere 
ozone is passed to the troposphere for eventual destruction with atmos- 
pheric constituents or the earth's surface (Craig, 1950, 1965). Craig 
further mentions that the process would have to be "most effective in 
middle and high latitudes in winter, when the total amount of ozone 
increases most rapidly" (Craig, 1950). 

Prabhakara (1963) performed a study of a mathematical model of 
ozone distribution based on photochemical equilibrium below 41 
kilometers and the effect of transport processes upon this distri- 
bution. His conclusion that "photochemical theory by itself cannot 
explain the latitudinal and seasonal variation of ozone" lends 
further support to ideas advanced by Craig (1950, 1965). 

In the past, sufficient ozone measurements on a hemispheric and 
vertical scale were unavailable. This prevented direct calculations 
of the transport of ozone by dynamic processes in the stratosphere. 

Newell (1961, 1962, 1964) advanced the idea that since the total ozone 
measured at the surface was highly correlated with the ozone amount in 
the 12 to 24 kilometer layer, then a measurement of the flux of total 
ozone would be indicative of the flux in the lower stratosphere. Using 
100 and 50 mb winds and total ozone measurements during the Inter- 
national Geophysical Year, Newell computed the mean and eddy fluxes 
of ozone in this layer. Newell concluded that the mean meridional 
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circulation and standing eddies did not serve as effective transport 
mechanisms. Thus, he argued that most ozone transport was attributed 
to quas i -horizonta 1 transient eddies. 

To explain the observed change in total ozone measured at the surface, 
a mechanism was proposed by Reed (1950) and Normand (1953). Although 
proposed separately, this theory is usually called the "Reed-Normand 
effect". The mechanism combines on a synoptic scale the horizontal 
advection and vertical motion to explain the variation of total ozone 
with the synoptic features of the upper troposphere and lower stratos- 
phere. Basically the argument is as follows. 

Air blows through the ridge and trough patterns of the height field. 
The temperature pattern associated with the height field in the lower 
stratosphere is fairly constant with warm areas in the base of troughs 
and cold areas in the ridges. Therefore, the quas i -s ta tionary tempera- 
ture patterns must be maintained by adiabatic warming and cooling of 
air parcels moving through them. The subsidence associated with the 
warming brings down ozone -enriched air from above. The lifting, 
associated with the cooling, transports ozone-poor air from lower 
altitudes. Since this region is below the level of the observed 
equatorial maximum, higher ozone values are brought southward to the 
trough and lower values are brought northward from equatorial regions 
to the ridge. 

The basic qualitative theories of movement and approximate distri- 
bution of ozone in the lower stratosphere have been proposed. Quanti- 
tative conformation of these theories has been lacking due to a 
minimum of information with respect to the horizontal and vertical 
distribution of ozone in the lower stratosphere (and elsewhere). With 
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the challenge thus defined, this study was undertaken to compare 
earlier theories with a quantitative description of ozone movements. 

As a first step the distribution of ozone mixing ratio on standard 
pressure surfaces in the lower stratosphere was constructed from the 
observations taken by the Air Force Cambridge Research Laboratories 
ozonesonde network during 29 April 1963 to 10 May 1963. It follows 
that ozone motions can be evaluated using basic meteorological data. 

To accomplish this evaluation, a computational scheme was proposed 
by Professor J. D. Mahlman of the Department of Meteorology, Naval 
Postgraduate School. This scheme puts the continuity equation for 
ozone (Haltiner and Martin, 1937) into a form that enables com- 
putation of the mean and eddy flux of ozone mixing ratio into a volume 
bounded by two pressure surfaces. The intent of the study is to perform 
quantitative calculations of ozone transport due to various physical 
processes based upon this formulation and assess their significance in 
terms of previous estimates. 
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II. PREPARATION OF DATA AND THE COMPUTATIONAL SCHEME 



A. BACKGROUND 

In 1963 an attempt was made to gather a more complete representation 
of the distribution of ozone on a synoptic scale. The Air Force 
Cambridge Research Laboratories (AFCRL) operated an eleven-station 
ozonesonde network in North America. Daily observations obtained by 
this network between 29 April 1963 and 10 May 1963 provide the basic 
data for this study (Hering and Borden, 1964). This period is just 
after the observed spring maximum of total ozone. Wind and temperature 
fields at standard levels were available from the National Weather 
Records Center, Asheville, North Carolina and from the Free University 
of Berlin (1963, 1966). 

B. PHYSICAL LIMITATIONS 

Considering the geographical location of the ozonesonde stations 
on the North American continent, an upper boundary of 24 kilometers 
defined by sufficient ozonesonde observations, and a lower boundary 
defined by the tropopause, the volume in which calculations can be 
made is limited. 

Four pressure levels - 200, 100, 50 and 30 mb - were chosen on 
the basis of readily available wind and temperature data. The 
selection of the grid boundaries at 30N-140W, 75N-140W, 75N-60W, and 
30N-60W minimizes the Atlantic Ocean, Pacific Ocean, and Northern 
Canada data -sparse areas. 

A la t i tilde- longi tude oriented grid was used in order to measure 
fluxes across boundaries of meridians and parallels. The grid 
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interval on a pressure surface is five degrees, which is Sufficient to 
represent synoptic scale motions in the stratosphere. 

For the four isobaric surfaces involved the distance between them 
is on the order of four kilometers. Considering the vertical distri- 
bution of ozone this interval is larger than desired and forces rough 
approximations in computing vertical derivatives. 

The time interval was chosen as one day corresponding to the 
frequency of ozonesonde observations. Considering the expected mcve- 
ment of ozone in this region and the synoptic scale of the analysis, 
this interval is reasonable. 

C. BASIC ATMOSPHERIC PARAMETERS 

Hand analyses of streamline, isotach, and isotherm patterns were 
prepared at each pressure surface for the time period of the calcu- 
lations. The daily 200 and 100 mb 1200 GMT Weather Bureau synoptic 
analyses and the daily 50 and 30 mb 0000 GMT synoptic analyses from 
the Free University of Berlin were used as a plotted data source. 

As shown in Appendix A, the streamline analyses describe typical 
patterns for this time of year. The lower levels reflect deep 
tropospheric systems that weaken with altitude. At the 30 mb surface 
the weakening winter circumpolar vortex is observed with weak anti- 
cyclonic centers forming and dissipating at the southern latitudes. 

The ten-day period of this investigation's major calculation is 
during the spring reversal of the upper stratospheric circulation 
from winter westerlies to summer easterlies* 

Since the ozonesonde observations were taken at 1200 GMT, an 
approximation to the parameters of the atmosphere at the 50 and 30 



21 



mb surfaces was made as follows. The zonal and meridional components 
of the wind field and the temperature field were linearly interpolated 
for 1200 GMT values from the 0000 GMT fields. 

Vertical motions were computed by solving the first law of 
thermodynamics in the form 
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The diabatic heating term was evaluated using the net heating rates 
as proposed by Kennedy (1964) for the spring period. Unrepresentative 
changes in vertical motions between grid points and excessive magni- 
tudes at certain grid points are the result of computational 
inaccuracies attributed to the assumption of a linear change of wind 
and temperature across a horizontal grid interval. To avoid these 
unrepresentative and dynamically untenable vertical motions, a 
smoothing scheme was used as follows : 

03. . = .15cjd. . + .15cjo. . + . 15o). , . +.15 cd._ l1 . + .4o>. .. (2) 

i,j i,j+i i>j-i i-i>j i+i, j i,j 

To aid in establishing the state of the atmosphere for the area 
of interest during the period of the calculations, ten day latitudinal 
averages, by pressure surface, were prepared for the temperature field 
and the zonal, meridional, and vertical components of the wind fields. 
These values are presented in meridional cross section format as 
figures 1-4. 

The values in figures 1-4 are in general agreement with results 
presented by other authors (e.g. Craig, 1965). The temperature field 
(fig. 1) is typical for the lower stratosphere at this time of year. 
Cold temperatures at southern latitudes give wav to the warm 
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LATITUDE 

FIG. 1. Ten-day Zonal Mean Temperatures (deg K) for the period 30 April to 
9 May 1963. Computed between longitude 60W and 140W. 
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LATITUDE 

FIG. 2. Ten-day Zonal Mean u-component of Wind (Knots) for the Period 
30 April to 9 May 1963. Computed between longitude 60W and 140W. 



o 



CN 




25 



LATITUDE 

FIG, 3. Ten-day Zonal Mean v-component of Wind (Knots) for the Period 
30 April to 9 May 1963. Computed between longitudes 60W and 140W. 
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temperatures associated with the middle and northern latitude troughs. 
This latitudinal temperature gradient weakens with altitude. The 
u-component of the wind (fig. 2) supports the earlier observation that 
the winter circumpolar vortex is weakening and being replaced from 
above by summer easterlies as can be seen at 30 mb and 35N. The mid- 
latitude wind maximum associated with the mean polar front jet position 
is seen at 200 mb and 45N. The v-component of the wind (fig. 3) shows 
little motion at the 30 and 50 mb levels. At 100 mb strong southerly 
movement north of 65N and strong northerly movement south of 50N is 
observed. At 200 mb, strong southerly motions are seen north of 65N. 

The vertical motions (fig. 4) show a general subsidence over the volume 
of calculations except for a weak rising area at 50 mb and strong upward 
motions between 45N and 50N at 200 mb. 

Most investigators presenting similar results use time scales of 
months or seasons. This study is limited to ten days and hence the 
values are not directly comparable. In addition, to properly discuss 
the various details of these mean atmospheric parameters, a separate 
and more thorough treatment is required than can be presented here. 
Therefore, further discussion of figures 1 - 4 is not presented. 

D. OZONE DATA AND ANALYSIS 

1 . The AFCRL Ozones onde Network 

For most of its duration the AFCRL ozonesonde network made 
weekly observations on Wednesday except for selected periods when more 
frequent observations were made. During one such period from 29 April 
1963 to 10 May 1963, daily observations were taken at 1200 GMT. 
Participating stations are listed in Table I. 
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OZONES ONDE NETWORK 



(after Hering and Borden, 1964) 



STATION 


LATITUDE 
(Deg. N) 


LONGITUDE 
(Deg. W) 


Albrook Field, Canal Zone (AWS) 


9.0 


79.6 


Colorado State University, Fort Collins 


40.6 


105.1 


Eielson AFB, Fairbanks, Alaska (AWS) 


64.8 


147.9 


Florida State University, Tallahassee 


30.4 


84.3 


Fort Churchill, Manitoba (Canadian Met. Br . ) 


58.8 


94.1 


Goose Bay, Labrador (Canadian Met. Br.) 


53.3 


60.4 


L. G. Hanscom Field, Bedford, Massachusetts 


42.5 


71.3 


Thule AFB, Greenland (AWS) 


76.5 


68.8 


University of New Mexico, Albuquerque 


35.0 


106.6 


University of Washington, Seattle 


47.4 


122.3 


University of Wisconsin, Madison 


43.1 


89.4 



TABLE I 
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The instrument used in determining the amount of ozone in the 
atmosphere was the dry chemiluminescent ozonesonde developed by Regener 
(I960), The ozonesonde was attached to a sounding balloon with a 
standard radiosonde package. The instrument measures the amount of 
ozone by drawing ozone-laden air across a disk coated with a chemical 
that produces a luminescent reaction with ozone. The intensity of the 
reaction is measured by a photomultiplier tube and in turn is tele- 
metered back to ground using standard radiosonde equipment (Regener, 
1960, 1964). Samples were taken every 15 seconds for a vertical 
resolution of approximately 250 feet (Hering and Borden, 1964). 

2 . Reliability of Measured Chemiluminescent Ozonesonde Data 
In view of the uncertainties inherent with the measurement 
and calibration of the Regener ozonesonde it is worthwhile to note 
the major corrections applied to the soundings and to note the results 
of several intercomparison tests with other ozonesondes . 

The disk coated with the chemiluminescent substance is sensitive 
to light, moisture, and ozone. Therefore, as a check, calibration 
shortly before launch is required (Regener, 1964). Hering and Borden 
(1964) report that achievement of reliable preflight calibration was 
difficult. This in turn required a correction to the telemetered data. 
The correction was achieved by adjusting the integrated observed ozone 
density to a simultaneously observed total ozone measurement using a 
spectrophotometer. The ozonagrams published by Hering and Borden (1964) 
have this correction applied. 

Hering and Dutsch (1965) compared the results of simultaneous 
soundings of the Regener chemiluminescent ozonesonde and the bubbler- 
type electrochemical ozonesonde. Their conclusion is that "ozonesondes 
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of the types tested are capable of providing reliable high-resolution 
measurements of the vertical ozone distribution". 

More recently Komhyr, Grass, and Proulx (1968) made intercomparison 
tests in various combinations between the Regener chemiluminescent 
ozonesondes , carbon-iodine ozonesondes , and Brewer-Mast ozonesondes. 

The result of the comparisons was that, after corrections, the 
chemiluminescent ozonesonde compared with the other ozonesondes in 
the stratosphere and gave "relatively reliable" results. 

Hering and Borden (1964) considered the data obtained by the 
ozonesonde network as "provisional". Considering the results of the 
two intercomparison tests and the observation that other investigators 
(Berggren and Labitzke, 1966, 1968; Breiland, 1964, 1965, 1967, 1968; 
Craig, DeLuisi, and Sticksel, 1967; Hering, 1966; Hering, Touart , and 
Borden, 1967; and Sticksel, 1966) have published studies using the 
corrected ozonesonde data, the reliability of the ozonesonde soundings 
as published by Hering and Borden (1964) is accepted. 

3 . Preparation of Ozone Data for Analysis 

Although the ozone soundings were accepted there are several 
features of the soundings that had to be modified to be consistent 
with the scale of the intended calculations. On the ozonagram thin 
layers with large unrepresentative gradients and layers in which the 
ozone values fluctuated excessively and with apparent randomness were 
considered not representative of the synoptic scale distribution. 

These features of the vertical distribution may either be a failure 
of the instrument or a sma 1 le r-sca le ozone distribution which is 
presently below our ability to identify. To overcome these apparent 
inconsistencies a subjective smoothing of the individual soundings was 
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completed as the first step in analysis. This smoothing is not to 
be over-emphasized as extreme care was taken not to destroy the 
characteristics of the soundings. The severest smoothing was on the 
order of .3 kilometers on the height scale with 95 percent of the 
ozonagram values being unaffected by the smoothing process. 

With the smoothed values obtained, a time series for the twe lve - 
day period was prepared for each station. Linearly interpolated values 
were calculated for stations missing an observation but having an 
observation 24 hours prior to and after the missing observation. 

Linear interpolations were also made in the case of the Fairbanks 
station where observations were reported at 0000 GMT. 

To supplement the small number of observations for the area under 
consideration, hand trajectories were calculated for the 24-hour 
periods prior to and following an observation. This technique used 
12-hour streamline and isotach patterns to obtain the position of the 
parcel 24 hours before and after the observation time. Under conserva- 
tion of ozone mixing ratio (x) , these new points then possess a value 
of X corresponding to the value measured at the original point. An 
estimate of the vertical motion during the trajectory was made using 
an approximation derived from the Eulerian expansion of the total 
derivative of pressure with respect to time, the hydrostatic equation, 
and the first law of thermodynamics. With this estimate of the 
vertical motion applied to the initial ozonagram sounding, a corrected 
trajectory value was obtained. The trajectory calculations are justi- 
fiable assuming the conservative property of ozone mixing ratio. As 
linearly interpolated values were not used for trajectory computations, 
this calculation increased the number of data points between two and 
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three times the original number of ozone data points available for 
analysis on a pressure surface. 

4 . Hand Analysis Techniques 

The manner in which the data was to be analyzed presented some- 
what of a problem. Synoptic models of ozone distribution on a pressure 
surface have not been published to the author's knowledge--with the 
exception of the six analyses presented by Berggren and Labitzke (1968). 
Other investigators have published information about the synoptic 
distribution of ozone but none have presented distributions on a 
pressure surface. Pittock (1969) has investigated the synoptic 
climatology of ozone using statistical methods. His "preliminary 
results" using 546 Brewer-Mast ozonesonde soundings are primarily in 
tabular format classed in.o general synoptic conditions by season. 
Breiland (1964, 1965, 1967, 1968) studied the distribution of ozone 
on a sub-synoptic scale primarily using vertical cross sections, time 
sections, and individual AFCRL network soundings. Sticksel (1966) 
analyzed the vertical distribution of ozone with data collected at the 
Tallahassee station of the AFCRL network. Briggs and Roach (1963) 
and Danielsen (1968) have investigated the ozone distribution, along 
with other parameters, near the jet stream and the tropopause using 
aircraft observations. 

The first attempt to analyze the plotted data points was incon- 
clusive. Merely constructing isolines of ozone mixing ratio from 
plotted data points left too much to the imagination. This resulted 
in inconsistent analyses when compared to the work of the previously 
mentioned authors and which were also consistent in time. 



32 



To overcome this problem the appropriate streamline analysis was 
put under the plotted ozone data. With the streamlines to guide the 
analyst in sparse data areas the desired consistency was achieved. To 
elaborate further, the criteria for placement of isolines of ozone 
mixing ratio are listed below in order of priority: 

a) observed data at map time; 

b) data from trajectory calculations as corrected for vertical 
motion ; 

c) streamline pattern; 

d) vertical motion patterns. 

Justification for this priority is obvious in the case of the first 
two criteria. The choice of the streamline pattern as the third 
criteria followed by the vertical motion patterns is based on Reed's 
(1950) conclusion that the horizontal advection contributed 2/3 of 
the day-to-day change whereas the vertical motions contribute 1/3 of 
the daily change. 

5 . Ozone Distribution in the Lower Stratosphere 

Before discussing the individual ozone analyses it is worth- 
while to compare the twelve-day mean ozone mixing ratio distribution 
obtained here with distributions prepared by other investigators. 

The twelve-day mean ozone mixing ratio computed for the period 
29 April to 10 May 1967 is presented in figure 5. It is seen that 
maximum values at southern latitudes and high altitudes give way to 
lower values at northern latitudes. The latitudinal gradient reverses 
between the 30 and 50 mb surfaces with the lower altitudes having 
higher values at northern latitudes. The latitudinal gradient re- 
verses between the 30 and 50 mb surfaces with the lower altitudes having 
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LATITUDE 

FIG. 5. Twelve-day Zonal Mean Ozone Mixing Ratio (mg g' 1 ) for the period 
29 April to 10 May 1963. Computed between longitudes 60W and 140W. 



higher values poleward. The slope of the isolines between 50 and 200 
mb is markedly downward. In comparison the isolines between 50 and 30 
mb are quasi-horizontal between 35N and 50N and rise northward between 
50N and 70N. 

A direct comparison of ozone mixing ratio values may be made with 
the values presented by Hering (1966) for the spring of 1963-1964. 

The two figures are similar except slightly higher values are given 
here using the 29 April 1963 to 10 May 1963 data. Although Hering 
uses the same data source as this study, these are the only available, 
directly comparable, mean ozone mixing ratio values computed from 
vertical soundings. 

Bojkov (1969) presents an annual vertical ozone distribution as 
determined by the Umkehr method. The units of this distribution are 
micromillibars and not easily comparable in a quantitative sense. 
Qualitatively, the same conclusions revealed by the twelve-day mean 
values are observed in Bojkov's distribution. 

Taking into account the information presented by Bojkov and Hering, 
it is concluded that the twelve-day mean cross section is in agreement 
with their results and therefore lends credibility to the individual 
analyses from which the mean distribution was prepared. A representa- 
tive sample of the distribution of ozone mixing ratio on the 200, 100, 
50 and 30 mb pressure surfaces is presented in Appendix A. 

The patterns made by the 30 mb ozone mixing ratio isolines are 
dominated by the effects of the horizontal component of the wind. 
Ozone-rich air is associated with the northward flow of the sub- 
tropical anticyclones forming on the southern edge of a weakening 
winter circumpolar vortex. Ozone-poor air is associated with the 
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southward flow on the western side of the large scale troughs super- 
imposed on the circumpolar vortex. The two contrasting air masses 
meet in the middle latitudes and mix. The air, continuing to move 



eastward, then diverges with relatively ozone-poor air associated with 
the southward flow on the eastern side of the subtropical anticyclones. 
Correspondingly, relatively ozone-rich air is associated with the 
northward flow on the eastern side of the large scale troughs. 

If long meridional trajectories form in connection with these 
cyclonic and anticyclonic wind systems, unusually large values of ozone 
can be observed at northern latitude stations while lower values are 
observed several hundred miles east or west of that station. The 
converse holds true for southern latitude stations. The movement of 
these wind systems will affect the daily observations at a station in 
that large day-to-day fluctuations of ozone can be caused by the 
associated wind shift. 

By comparing the vertical motion fields on a pressure surface 
(not reproduced here) with the ozone mixing ratio distribution, some 
very interesting observations can be made. The western side of the 
anticyclones are associated with rising vertical motion, northerly 
trajectories, and high ozone values while the western side of troughs 
are associated with sinking motions, southerly trajectories, and low 
ozone values. Tims it is qualitatively observed that the Reed-Normand 
effect is operating in the opposite sense than at lower altitudes. 

The result of this mechanism is in evidence in the mean ozone mixing 
ratio distribution of figure 5 where the isolines between 30 and 30 
mb slope' upward at northern latitudes. 
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Masses of air with similar ozone characteristics can be identified 



and tracked for periods of several days. Such identifiable air masses 
seem to be prevalent in the wind fields of the southern anticyclones. 

Ozone mixing ratio patterns at 50 mb are similar to the patterns 
at 30 mb with several exceptions. Streamline patterns at this level 
reflect the upper, weakening portions of tropospheric systems. Weaker 
southern anticyclones are also present. The ozone mixing ratio pattern 
shows ozone-rich air associated with the southward flow that at 30 mb 
is associated with ozone-poor air. Similarly, ozone-poor air is 
associated with the northward flow that at 30 mb is associated with 
ozone-rich air. 

A qualitative comparison of the vertical motion patterns with the 
ozone mixing ratio patterns show the Reed-Normand effect to be opera- 
ting in the sense originally proposed by Reed (1950) and Normand (1953). 
The sinking, southward moving air into troughs are areas of high ozone 
values while rising northward moving air is associated with low ozone 
va lues . 

Masses of ozone-rich and ozone-poor air can be identified as at 
30 mb. However, at both pressure levels it noted that these identi- 
fiable ozone masses may be an unnatural consequence of the trajectory 
supplementation of the data and a weak latitudinal gradient of ozone. 

The ozone mixing ratio patterns at 100 and 200 mb are similar to 
the 50 mb patterns. The effect of the streamline pattern is more 
pronounced . 

The vertical motion pattern on a broad scale reflects the Reed- 
Normand effect; however, vertical motion maxima with horizontal 
dimensions on the order 2 or 3 grid intervals are not clearly associated 
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with the Reed-Normand predicted high and low ozone values. This is 
caused by the density of the ozone data being insufficient to isolate 
these smaller scale features or by bad a) calculations. Support for 
this point of view is gained from a calculation of the local change, 
horizontal advection and vertical advection terms of the Eulerian 
expansion of the tota 1 d eri va t ive of ozone mixing ratio with respect 
to time. This calculation was made on a point by point basis for the 
entire grid. The results were inconclusive since the error term, 
computed by assuming ozone is conserved and summing the three remain- 
ing terms, was too large in most cases. This indicated a smaller 
scale ozone distribution that was not able to be identified with the 
data available. 

The ozone analyses prepared for this study compare in a general 
manner with the distributions used by Breiland (1964, 1965, 1967, 
1968), Sticksel (1965), Pittock (1969), Briggs and Roach (1963), 
Danielsen (1968), and Berggren and Labitzke (1965, 1968) in their 
respective studies. Only the 200 mb analyses for 3 and 4 May 1963 
presented by Berggren and Labitzke (1968) are directly comparable 
with the analyses compiled here. The analyses are similar in that the 
major features are the same but minor differences occur with the 
smaller features. It is significant to note that although the AFCRL 
data was the initial input in both cases, Berggren and Labitzke used 
isentropic trajectories to supplement their data points while this 
study used 24-hour isobaric trajectories corrected for vertical 
mot i on . 

Willi the data apparently consistent with the earlier theories, 
previous investigators, and consistent with respect to other 
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atmospheric parameters it was decided that these analyses would provide 
an acceptable data base for a quantitative calculation of ozone trans- 
port mechanisms. 



E. THE OZONE BALANCE EQUATION FOR A VOLUME BOUNDED BY TWO PRESSURE 
SURFACES 



As a first step in evaluating earlier theories of ozone distri- 
bution and movement, the calculation of the mean and eddy fluxes of 
ozone on a synoptic scale are desirable. Considering the minimum 
amount of ozonesonde observations at a given time, the uncertainty in 
the synoptic distribution of ozone on a pressure surface, and the un- 
certainty in the performance of the ozonesonde, an averaging scheme 
over many' data points is preferred as compared to a point by point 
evaluation. To achieve this, the continuity equation for ozone is 
written in pressure coordinates as 
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(3) 



By putting the advection terms in equation (3) into flux form 
and taking the area average, equation (3) becomes 

f - - ^>4 - If + < P - *'> + i ( K p IP + (4) 

Now by applying the divergence theorem equation (4) may be 
rewritten as 
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(5) 



where c^, the component of the wind normal to the boundary, is defined 
as positive toward the center of the volume. 
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Along a line let X = X + X* where x is the line average and x' is 



the point deviation from the line average. Similarly for averaging 
operations, equation (5) becomes 
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noted 



it is assumed 
that fdl = A1 
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pressure surface and 



Now average equation (6) between two pressure surfaces, 
p , to obtain the ozone balance equation 



p and 
r u 



(MFLA) 


(EFLA) 


(MFLO) (MFT) 




A A VW- 

A1 Al — r 

r xc „ + r x c n - 


CD X 1 , CD X I 

Ap P L Ap p u 




(EFLO) 


(EFT) 


(PR) (SUBLO) 




0)*X* | . 0)*X* 1 

A P 'p L Ap p u 


Ap dp 'p L 


(7) 


(SUBT) 


(SUB LA) 







Sx | _ A l ( v dx\ 

Ap r ^P P u A. V H dn^’ 



where the vertical average has been applied to evaluate the terms 

in a form for numerical computations. 

The terms of the equation are : 

(LC) the time rate of change of average ozone mixing 

ratio over the volume; 

(MFLA) the mean flux of ozone into the region through the 

lateral boundaries (the vertical side of the volume); 

(EFLA) the eddy flux of ozone through the lateral boundaries 
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(MFLO) 


the mean flux of ozone into the region across 
bottom boundary; 


the 


(MFT) 


the mean flux of ozone into the region across 
top boundary; 


the 


(EFLO) 


the eddy flux of ozone into the region across 
bottom boundary; 


the 


(EFT) 


the eddy flux of ozone into the region across 
top boundary; 


the 


(PR) 


the average difference between production and 
recombination of ozone for the whole volume ; 




(SUBLO) 


the flux of ozone into the region, across the 
boundary, due to sub-grid scale diffusion; 


lower 


(SUBT) 


the flux of ozone into the region, across the 
boundary, due to sub-grid scale diffusion; 


upper 


(SUBIA) 


the flux of ozone into the region through the 
due to sub-grid scale diffusion. 


s ides 



All the terms have units of micrograms per gram per day. This 
unit is convenient in view of the time scale of the observations. In 
addition, the value of the ozone mixing ratio at a given height appears 
as a variable on the ozonagram. 

Computer programs were written to evaluate all terms on a daily 
basis for the 30/50, 50/100, and 100/200 mb volumes with the ex- 
ception of the production and recombination term and the sub-grid 
scale diffusion terms. 

Brewer and Wilson (1968) have evaluated the meridional net ozone 
production rate for the winter and summer season from the equator to 
sixty degrees latitude and between 200 and 30 mb. Since the evaluation 
of the PR term would require a complete and separate treatment. Brewer 
and Wilson's results are accepted for use in this calculation. To 
compute an estimated value for the net ozone production rate, a rough 
interpolation between seasons and pressure surfaces followed by an 
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estimated latitudinal average was necessary. The appropriate unit 
conversions were made. These values will be presented with the 
results of the ozone balance equation in a later section. 



The sub-grid scale diffusion terms are expected to be small when 
compared to the magnitude of the other terms. Lacking a reliable 
estimate of the eddy-diffusion coefficients of ozone mixing ratio, 
these terms were not calculated. 

The initial attempt to evaluate the ozone balance equation by 
using the observed wind to directly evaluate each term resulted in 
unusually large errors that could be traced to the fluctuations of 
the mean terms. To improve upon these results the mass continuity 
equation was solved for the mean normal component of the wind at the 
side boundaries and in terms of the computed cd fields. The evaluation 
of the mean normal component of the wind can be obtained by applying 
the divergence theorem to the area -averaged continuity equation and 
solving the following form, 



/VW' 

c = 
n 



A (cd - CD ) 



A1 ( P L ■ P u> ' 



( 8 ) 



Although this scheme assured the conservation of mass it increases 
the dependence of the calculations on the computed vertical motions 
as compared to directly observed winds. 
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III. EVALUATION OF THE TERMS OF THE OZONE BALANCE EQUATION 



Without further modification the terms of the ozone balance 
equation were computed. The daily calculations tended to fluctuate 
rapidly and may be unrealistic. Considering the magnitude of the error 
term in comparison with the magnitude of the other terms, it was decided 
to present the ten-day mean values of the terms. These values are given 
in tabular format in Table II, In the computed format, the individual 
mean and eddy terms cannot be compared directly. To effect a com- 
parison the mean terms must be summed thereby giving the net effect of 
the mean terms on the volume. 

To obtain a form in which the mean horizontal and vertical flux 
may be directly compared to the eddy flux terms, the sum of the mean 
terms, as presented in Table II, may be computed in another manner. 

The area-averaged continuity equation can be written as a repre- 
sentation of the net effect of the mean terms at the middle pressure 
between two pressure surfaces. By applying a vector identity and the 
divergence theorem, the first and third terms on the right hand side 
of equation (6) can be written as 




(9) 



By combining like terms and putting into computational format, 



equation (9) becomes 



(MLAD) 



(MVAD) 




( 10 ) 
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Key: See Section II-E 

SUM = MFLA + MFLO + MFT 

AS COMPUTED: 



Vo lume 


(LC) (MFLA) (EFLA) (MFLO) 


(MFT) 


(EFLO) (EFT) 


(PR) 


(ERROR) 


30/50 mb 


-173 = + 184 


+ 54 -688 


+696 


- 67 -50 


-125 


+177 


50/100 mb 


-517 = +1102 


- 46 -921 


+275 


-103 +27 


- 44 


+806 


100/200 mb 


-273 = - 229 


-321 - 47 


+461 


-155 +51 


- 16 


+ 17 


COMPARISON 


FORMAT: 












Volume 


(LC) (SUM) 


(EFLA) 


(EFLO) 


(EFT) 


(PR) 


(ERROR) 


30/50 mb 


-173 = +192 


+ 54 


- 67 


-50 


-125 


+ 177 


50/100 mb 


-517 = +456 


- 46 


-103 


+27 


- 44 


+806 


100/200 mb 


-273 = +185 


-321 


-155 


+51 


- 16 


+ 17 



TABLE II. The Ten-day Mean Ozone Balance for Indicated 
Volumes. Symbols for each term given in 
Section II-E. Units are micrograms per gram 
4 

per day x 10 . 
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The two terms are defined as: 



(MLAD) the mean advection through the lateral boundary; 

(MVAD) the advection into the volume due to the mean 

vertical motion. 

The two terms of equation (10) were evaluated using the same input 
values as for the evaluation of equation (7). The results of the 
evaluation of equation (10) are presented in pictorial form in figure 
6 with the local change and eddy terms as evaluated by equation (7). 
The net production term is the same in both cases. The error terms 
in each evaluation are different only to a small degree. It is 
therefore concluded that the computational scheme used to arrive at 
the value of the mean terms is a reasonable evaluation of the input 
parameters . 

With the mean and eddy ozone flux terms in a format for comparison 
and some degree of confidence in the method of evaluation of the mean 
terms, it is appropriate to discuss the individual terms of the ozone 
balance as presented in figure 6. 

A. THE LOCAL CHANGE TERM 

Computation of this term was made directly from the distribution 
of ozone on a pressure surface and therefore is not affected by 
other computed parameters. A rough check on the value of this term 
was made by comparing the computed result, multiplied by 100, with a 
100-day difference between area averages estimated from Godson's 
(1960) mean distribution of total ozone. This period falls after 
the March maximum of total ozone shown by Godson and the sign of the 
mean local change term is negative. With these observations in mind, 
the 100-day period between May and August was chosen for comparison. 
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FIG. 6. Ten-day mean ozone balance in 
advective format. Symbols for each term 
given in Section II-E and Section III. ^ 
Units are micrograms per gram per day x 10 . 
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After making the appropriate unit conversions, the mean local 
change term accounted for 85 percent of the change in Godson's distri- 
bution. Thus the rate of decrease of ozone within the volume for this 
period appears representative of the seasonal trend. 

It is significant to note that the maximum of total ozone in the 
atmosphere is observed in March and a net hemispheric decrease is ob- 
served after that month (Godson, 1960). The dynamic mechanisms or 
conditions that created the maximum conditions, particularly at the 
northern latitudes, must have ceased operating. One may speculate 
that the mechanism to accomplish this winter-spring northern maximum 
is the long meridional trajectories set up by the breakdown of the 
polar night vortex associated with the mid-winter "sudden warmings". 
Further, the dynamic meachnisms or conditions now operating must be 
serving to deplete the atmosphere of ozone. 

B. THE MEAN TERMS 

The mean advection through the lateral boundaries can only be 
evaluated as to direction into or out of the volume. In view of this 
limitation, the only point to be made is a mention that the term in 
each volume is contributing a small net loss to the volume and the 
mean terms (MLAD) are smaller than their respective eddy terms (EFLA) 
(see fig. 6). 

On the other hand, the mean advection in the vertical shows a 
consistent downward motion through the three volumes of interest. In 
the two upper volumes the mean vertical advection is the dominant 
term. 
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c. the eddy terms 



In the 30/50 mb volume the lateral eddy flux term (EFLA) is an 
input to the volume and agrees with the earlier qualitative observation 
that the anticyclones on the southern edge of the circumpolar vortex 
appeared to be transferring ozone northward. The magnitude of the 
lateral eddy flux term (EFLA) is small relative to the observed change 
within the volume. 

The lateral eddy flux terms at 50/100 and 100/200 mb serve to 
deplete the volumes of ozone. The value of the term at 50/100 mb is 
small relative to the mean vertical advection and the vertical eddy 
flux across the 100 mb surface. In the 100/200 mb volume the lateral 
eddy flux term is the largest term in the ozone balance. An interest- 
ing and unexplained question is the ultimate destination of the ozone 
advected out of the volume of computations by the lateral eddy flux. 

The vertical eddy flux terms show a consistent downward motion, 
from 50 to 200 mb and are a measure of the vertical motion portion 
of t lie Reed-Normand effect. At these levels the direction is con- 
sistent with that proposed by Reed and Normand. At the 30 mb surface 
the direction of the vertical eddy flux is upwards. This can be 
traced to the northward, rising ozone-rich air associated with the 
western side of the anticyclones at this level. This is consistent 
with the earlier section in which the Reed-Normand effect was observed 
to be operating in the opposite sense as observed in lower altitudes. 

By comparing the vertical eddy flux at the upper and lower sur- 
faces at each volume, it can be seen that these eddies contribute a 
net loss of ozone to each volume. Another interesting, unanswered 
question is posed by this observation. If a given amount of ozone is 
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going downward into a trough, what physical process or dynamic mechanisms 
are operating to prevent an equal amount of ozone from proceeding back 
upwards to the ridge? The answer probably lies with a smaller scale 
turbulence phenomenon not yet identified much as the energy cascade 
in the lowest level of the atmosphere. 

The vertical eddy flux terms are the most significant of all the 
terms computed in the ozone balance. This conclusion is attributed 
to the fact that the vertical eddy flux computed for this particular 
area's pressure surfaces may be representative of the whole global 
pressure surface in question. This cannot be said of the lateral 
eddy flux terms. As is to be described in a later section, the 
magnitude of the mean terms is subject to question. 

D. THE PRODUCTION AND RECOMBINATION TERM 

In comparison with the magnitudes of the other terms, the pro- 
duction and recombination term as estimated from Brewer and Wilson*s 
(1968) values does not have a significant effect except in the 30/50 
mb volume. Here the magnitude is about 2/3 of the local rate of 
change. In future investigations this must be taken into con- 
sideration . 

E. THE ERROR TERM 

The error term in Table II and figure 6 is the value that must be 
added to the local change in order to equal the sum of the mean, eddy 
and net production terms. In numerical computations of this sort, 
errors are expected because of data uncertainties and since many 
approximations are made to facilitate the evaluations of the individual 
terms and their components. 
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The mean errors in each volume are of the order of magnitude of the 
largest term. The mean errors at 30/50 and 50/100 mb are representative 
of the daily errors at these levels from which the mean error was derived. 
This is not true at the 100/200 mb volume where the daily errors average 
one order of magnitude larger than the ten-day mean error. Thus it is 
noted that the ten-day mean error is unrepresentative of the error in 
the 100/200 mb volume and, further, that the calculations here suffer 
from the same problems as the upper two volumes. 

Considering that the error in all three volumes indicates too much 
ozone on the right hand side of the balance, the failure to compute the 
sub-grid scale diffusion terms has to be considered. It is expected 
that these terms will evacuate ozone from the volume with a magnitude 
considerably smaller than the magnitude of the error. Therefore it 
is improbable that the failure to calculate these terms contributes 
much to the error. 

The first serious possible source of error is that the ozone 
mixing ratio distributions on the pressure surfaces are not accurate 
enough for the calculations. Only by having a more dense data network 
will this ever be known for sure. The probability that this is the 
major source of error is diminished considering the consistency in 
daily analyses and the consistency of the distributions with the 
theories and distributions advanced by earlier investigators. 

The second and most probable source of error is in the evaluation 
of the vertical motions. As was noted earlier the evaluation of the 
ozone balance equation is heavily dependent upon vertical motions. 

While a consistent error in the vertical motions will not affect the 
eddy terms greatly, it will have considerable effect on the mean terms. 
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Since the continuity equation was used to evaluate the mean flux 
through the lateral boundaries, this term is affected also. This line 
of reasoning is consistent with the fact that the mean terms in most 
cases are large compared to the rest of the terms as can be seen in 
Table II and figure 6. It can be reasoned that a consistent error in 
the vertical motion at a particular pressure level could affect each of 
the mean terms enough to produce a significant error. 

The size of the error in the vertical motion does not have to be 
very large to seriously affect the ca lculations . To see this, consider 
the vertical motion scheme used in the computations and presented as 
equation (1). The diabatic heating term is the term most subject to 
question as its values determined by any method are highly uncertain. 
This is attributed to the lack of an accurate method to check the 
results of the determination of the diabatic heating. If the area 
average of equation (1) is taken and it is assumed that the cor- 
rection to the vertical motion (a> ) is a function only of the 

correction to the diabatic heating term (corr) then it can be 
written that approximately, 



0) 

corr 



-corr 

a/c 

P 



(ID 



A representative value of area averaged diabatic heating at 
75 mb as estimated from Kennedy's (1965) results can be considered 
-.35 degrees per day. If this were incorrect by an arbitrarily 
selected + .2 degrees per day, then the evaluation of equation (11) 
using representative values for evaluation of the specific volume 
would produce 

^corr = — *24 mb per day, or about a 25 percent error in the mean 
vertical motions. 
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The mean vertical advection of ozone across a pressure surface is 

represented by cjd in equation (10). Evaluating the effect of a) 

op corr 

on this term using representative values of X produces uncertainty in 

-4 

this term of approximately 210 x 10 micrograms per gram per day. It 
is easily seen from equations (7) and (8) how this would affect the 
mean flux through the lateral boundaries as presented in Table II. 

The value of the mean advection through the lateral boundary is small 
(fig. 6). Therefore, it is reasoned that in this form the mean lateral 
advection is not a major source of error. Similar arguments can be 
made for the 30/50 and 100/200 mb mean terms. 

It can be concluded that the uncertainty in the evaluation of 
mean stratospheric vertical motions is large enough to result in the 
observed magnitude of the error to the ozone balance equation through 
evaluation of the mean terms. Implied from the excess of o zone 
represented by the error term is that the net input to the particular 
volume by the mean terms, particularly the mean vertical flux (or 
advection) terms, is too large although the sense of the individual 
mean terms is probably correct. 
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IV. ADDITIONAL CALCULATIONS 



An attempt was made to support the results of the qualitative 
evaluation of the distributions of ozone on a pressure surface and 
the terms of the ozone balance by calculating the correlation 
coefficients of the meridional wind and ozone mixing ratio and the 
correlation coefficients of the vertical motion and the ozone mixing 
ratio. The calculations are presented as figures 7 and 8. For the 
number of data points involved a correlation coefficient greater 
than .2 is statistically significant at the 95 percent level of 
c onf idence . 

The correlation coefficient for two variables is defined as the 
covariance divided by the variance. It follows then that the sign 
of the correlation coefficient is also the sign of the covariance. 

This relationship allows comparison with the eddy terms of the ozone 
balance . 

Looking at figure 7, it can be seen that there are non-significant 
values of the correlation coefficients for most of the 30 and 50 mb 
surfaces during the ten-day period. This is consistent with the small 
values of the lateral transport computed by the ozone balance. The 
ozone distributions observed on the pressure surfaces are evaluated 
as having been spread out by the weak wind fields at these levels and 
accomplishing little transport. The positive area at 30 mb and 70N 
latitude is not conclusive since this area was determined by one data 
point. It is consistent, however, with the assumption that ozone-poor 
air is advected southward by the horizontal eddy flow at this altitude. 
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LATITUDE 

FIG. 7* Ten-day Mean Correlation Coefficients Between Meridional Wind 
and Ozone Mixing Ratio for the Period 30 April to 9 May 1963. Computed 
between longitude 60W and 140W. 



At the 100 mb level the negative values to the north are associated 
with high ozone values being moved southward by the horizontal eddy 
flow. The positive values to the south are an indication of a counter- 
gradient transport (see fig. 5) of high ozone values northward and low 
ozone values southward by the horizontal eddy flow. To evaluate the 
cause of the counter-gradient transport, the vertical motion patterns 
(not presented here) were compared with the streamline patterns. It 
was observed that in many cases the maximum sinking areas connected 
with ozone-rich air were slightly (45 degrees in phase angle) to the 
east of the base of the trough. This offset of the vertical motion 
pattern from the trough line has been discussed by Ohring and Muench 
(1960). The sinking consequently offset the entire ozone mixing ratio 
pattern to the east and set the conditions for counter-gradient flow 
by the horizontal eddies. 

The results at 200 mb were expected to be more conclusive than 
shown in figure 7. The values reflect the stronger relationships at 
100 mb, however, they are not at the 95 percent level of significance. 
Without further calculations, a definite conclusion cannot be advanced. 

The results of figure 8 strongly support earlier contentions. The 
weak positive values are associated with upward eddy transport of high 
ozone values by the anticyclonic activity at 30 and 50 mbs. The high 
negative values at the southern latitudes indicate the expected down- 
ward transport of high ozone values. These latitudes coincide with 
the base of the troughs on the circumpolar vortex and hence general 
areas of maximum sinking motions. At 100 mb the secondary maximum of 
significant negative values at the northern latitudes are evaluated as 
downward transport of high ozone values by the weaker semi-permanent 
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trough of the circumpolar vortex. Based on the values at 200 mb it is 
further reasoned that the vertical eddy flux of ozone across the 
tropopause is taking place in connection with the maximum sinking 
areas near the base of the transient troughs on the circumpolar 
vortex . 

In the recent past several investigators (Newell 1961, 1962, 1963, 
1964; Mahlman, 1966, 1967; and Molla and Loisel, 1962) have used the 
latitudinal covariance of the meridional and vertical wind components 
(v'w 1 ) as an indication of stratospheric eddy transport processes. As 
a further check on the validity of the terms of the ozone balance 
equation, it was decided to compare (v'w*), by using correlation 
coefficients, with the information prepared for this study. This 
distribution is presented in figure 9. 

The negative values of the correlation coefficients are associated 
with a northward and downward transport while the positive areas 
indicate a southward and downward transport. The information from 
fig. 9 is in agreement with the sense of the vertical eddy flux terms 
of the ozone balance, and the ozone distribution on a pressure surface 
except at the 30 mb surface. At this level the mean ozone mixing 
ratio, streamline and vertical motion patterns show a northward and 
upward flow of high ozone values. This was confirmed by a vertical 
eddy flux out of the 30/50 mb volume across the 30 mb surface and 
figure 8. The vertical eddy flux is a direct calculation of transport 
as compared to the correlation coefficients which are only an indi- 
cation of transport. This implies that the correlation coefficient 
must be on the order of magnitude of .5 or greater to give confidence 
in the conclusion that this term is in fact an indicator of transport. 
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FIG* 9. Ten-day Mean Correlation Coefficients Between Meridional Wind 
and Vertical Motion for the Period 30 April to 9 May 1963. Computed 
between longitude 60W and 140W. 



V. SUMMARY AND CONCLUSIONS 



Even though the magnitude of the mean terms is subject to question, 
the distribution of ozone on a pressure surface, the ozone balance, and 
the supporting calculations allow several general conclusions to be 
drawn. 

1. It is possible to construct the distribution of ozone on a 
pressure surface from the Regener ozonesonde data and then perform 
numerical calculations that provide reasonable results. 

2. During this time period, the atmosphere is past its yearly 
maximum content of ozone at northern latitudes. There is no evidence 
to support a large northward movement of ozone at this time. 

3. During the ten-day period in question, ozone content of the 
volume is decreasing. To accomplish this, the mean vertical motions, 
supported by the vertical eddy transport at 50, 100, and 200 mb, are 
moving ozone downward to the tropopause level where the vertical eddy 
transport associated with the transient eddies on the circumpolar 
vortex then move the ozone into the troposphere for eventual 
destruction. 

4. The vertical eddy transport is a measure of the vertical part 
of the Reed-Normand effect. This mechanism is operating in the sense 
as proposed between 50 and 200 mb. At 30 mb the Reed-Normand effect 
is operating in the opposite sense with ozone-rich air being advected 
upward . 
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APPENDIX A 



STREAMLINE AND OZONE MIXING RATIO ANALYSIS FOR 3 MAY 1967 

This appendix presents a representative sample of the distributions 
of ozone mixing ratio on a pressure surface in the lower stratosphere. 
To understand the effect of the wind field upon the ozone distribution, 
the appropriate streamline pattern is presented. 

The choice of date was made to facilitate comparison with the pat- 
terns published by Berggren and Labitzke (1968). 
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FIG. A-l. The 200 mb Streamline Analysis for 3 May 1963. 
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FIG. A-3. The 100 mb Streamline Analysis for 3 May 1963. 
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The 100 mb Ozone Mixing Ratio Analysis (mg g 
for 3 May 1963. 
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FIG. A-5. The 50 mb Streamline Analysis for 3 May 1963. 
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FIG. A-6. The 50 mb Ozone Mixing Rat 
for 3 May 1963. 
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FIG. A-7 . The 30 mb Streamline Analysis for 3 May 1963. 
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FIG 0 A-8. The 30 mb Ozone Mixing Ra 
for 3 May 1963. 



LIST OF REFERENCES 



Berggren, R. and K. Labitzke, "A detailed study of the horizontal and 
vertical distribution of ozone, " Te llus , v. 18, p. 761-722, 1966. 

, "The distribution of ozone on pressure surfaces," Te llus , 

v. 20, p. 88-97, 1968. 

Bojkov, R. D., "Vertical distribution of ozone in the earth's atmos- 
phere," Meteorologua i Gedrologua , v. 10, p. 3-11, 1955. 

(Translated by AFCRL, G.T.R. 66-22, Oct 1966.) 

, "Computing the vertical ozone distribution from its relation- 
ship with total ozone amount," Journal of Applied Meteorology , 
v. 8, p. 284-292, April 1969. 

Breiland; J. C., "Vertical distribution of atmospheric ozone and its 
relation to synoptic meteorological conditions," Journal of 
Geophysical Research , v. 69, p. 3801-3808, 15 September 1964. 

, "A case study of the vertical distribution of atmospheric 

ozone," Journal of Applied Meteorology , v. 4, p. 357-364, June 
1965. 

, "Compar is.on of the vertical distribution of thermal stability 

in the lower stratosphere with the vertical distribution of 
atmospheric ozone, 1 ' Journal of the Atmospheric Sciences , v. 24, 
p. 569-576, September 1967. 

y "Some large scale features of the vertical distribution of 

atmospheric ozone associated with the thermal structure of the 
atmosphere," Journal of Geophysical Research , v. 73, p. 5021-5028, 
15 August 1968. 

Brewer, A. W.* and A. W. Wilson, "The regions of formation of atmospheric 
ozone," Quarterly Journal of the Royal Meteorological Society , 
v. 94, p. 249-265, July 1968. 

Briggs, J. and W. T. Roach, "Aircraft observations near jet streams," 
Quarterly Journal of the Royal Meteorological Society , v. 89, 
p. 225-247, April 1963. 

Craig, R. A., The Observations and Photochemistry of Atmospheric Ozone 
and Their Meteorological Significance , American Meteorological 
Society, September 1950. 

, The Upper Atmosphere: Meteorology and Physics , Academic 

Press, 1965. 



69 



, J. D. DeLuisi, and P. R. Sticksel, "Ozone distribution over 

Tallahassee, Florida," Journal of Geophysical Research , v. 72, 
p. 1661-1665, 15 March 1967. 

Danielsen, E. F., "Stratospheric- tropos pheric exchange based on radio- 
activity, ozone and potential vorticity," Journal of the Atmos- 
pheric Sciences , v. 35, p. 502-518, May 1968. 

Dobson, G. M. B. and D. Harrison, "Measurements of the amount of ozone 
in the earth’s atmosphere and its relation to other geophysical 
conditions," Proceedings of the Royal Society of London , v. 110A, 
p. 660-693, April 1926. 

, D. Harrison, and J. Lawrence, "Measurements of the amount of 

ozone in the earth’s atmosphere and its relation to other geo- 
physical conditions - Part II," Proceedings of the Royal Society 
of London , v. 114A, p. 521-541, 1927. 

, D. Harrison, and J. Lawrence, ’'Measurements of the amount of 

ozone in the earth’s atmosphere and its relation to other geo- 
physical conditions - Part III," Proceedings of the Royal Society 
of London , v. 122A, p. 456-486, February 1929. 

, "Observations of the amount of ozone in the earth’s atmosphere 

and its relation to other geophysical conditions - Part IV," 
Proceedings of the Royal Society of London , v. 129A, p. 411-433, 
November 1930. 

Free University of Berlin, Tagliche Hohenkarten der 50-mb-Flache sowie 
monatliche Mittelkarten fur das Jahr 1963 , by Ute Hengelhaupt, 
Ingrid Reuter, Richard Scherhag, Gunter Warnecke and others, v. 38, 
p. IV/29-V/10, 1966. 

, Daily and Monthly Northern Hemisphere 30-Millibar Synoptic 

Weather Maps of the Year 1963 , by Barbara Kriester, Karin 
Labitzke, Drazen Poje , Richard Scherhag, and Renate Stuhrmann, 
v. 39, p. IV/29-V/10, 1963. 

Godson, W. L. , "Total ozone and the middle stratosphere over arctic and 
sub-arctic areas in winter and spring," Quarterly Journal of the 
Royal Meteorological Society , v. 86, p. 301-317, July 1960. 

Gotz , F. W. P. , "Ozone in the atmosphere," Compendium of Meteorology , 
p. 275-291, Waverly Press, 1951. 

Haltiner, G. J. and F. L. Martin, Dynamical and Physical Meteorology , 
McGraw-Hill Book Company, 1957. 

llering, W. S., "Ozone and atmospheric transport processes," Te llus , 
v. 18, p. 329-336, 1966. 



70 



, and T. R. Borden Jr., Ozonesonde Observations Over North America , 

Volume 2 , Environmental Research Papers, No. 38, Air Force Cambridge 
Research Laboratories, July 1964. 

, and T. R. Borden Jr., Mean Distributions of Ozone Density Over 

North America, 1963-1964 , Environmental Research Papers, No. 162, 

Air Force Cambridge Research Laboratories, December 1965. 

, and H. Dutsch, "Comparison of chemiluminescent and electro- 
chemical ozonesonde observations , 11 Journal of Geophysical Research , 
v. 70, p. 5483-5490, 15 November 1965. 

, C. N. Touart, and T. R. Borden Jr., "Ozone heating and radiative 

equilibrium in the lower stratosphere ," Journal of the Atmospheric 
Sciences , v. 24, p. 402-413, July 1967. 

Kennedy, J. S., Energy Generation Through Radiative Processes in the 

Lower Stratosphere , Planetary Circulations Project, Report No. 11, 
Massachusetts Institute of Technology, 1964. 

Komhyr, W. D. , R. D. Grass, and R. A. Proulx, Ozonesonde Intercomparis on 
Tests , ESSA Technical Report ERL 85-APCL4, September 1968. 

Kulkarni, R. N. , "Comparison of ozone variations and of its distri- 
bution with height over middle latitudes of the two hemispheres," 
Quarterly Journal of the Royal Meteorological Society , v. 88, 
p. 522-534, October 1962. 

Mahlman, J. D., The Atmospheric General Circulation and Transport of 
Radioactive Debris , Atmospheric Science Paper No. 103, Colorado 
State University, 1966. 

, Further Studies on Atmospheric General Circulation and Trans- 
port of Radioactive Debris , Atmospheric Science Paper No. 110, 
Colorado State University, 1967. 

Mateer, C. L. and W. L. Godson, "The vertical distribution of atmos- 
pheric ozone over Canadian stations from umkehr observations," 
Quarterly Journal of the Royal Meteorological Society , v. 86, 
p. 512-518, October 1960. 

Meetham, A. R. , "The correlation of the amount of ozone with other 
characteristics of the atmosphere," Quarterly Journal of the 
Royal Meteorological Society , v. 63, p. 289-307, 1937. 

Molla, A. C. and C. J. Loisel, "On the hemispheric correlations of 
vertical and meridional wind components," Geofisica Pura E 
Applicata , v. 51, p. 166-170, 1962. 

Newell, R. E., "The transport of trace substances in the atmosphere and 
the implications for the general circulation of the stratosphere," 
Geofisica Pura E Applicata , v. 49, p. 137-158, 1961 



71 



, n Transfer through the Tropopause and Within the Stra t os phere , " 

Planetary Circulations Project, Report No. 7, Massachusetts 
Institute of Technology, 25 October 1962. 

, ’’The general circulation of the atn^osphere and its effect on 

the movement of trace substances,*' Journal of Geophysical Research , 
v. 68, p. 3949-3962, 1963. 

, "Further ozone transport calculations and the spring maximum 

in ozone amount," Geophysics, Pure and Applied , v. 59, p. 191-206, 
1964. 

Normand, Sir Charles, "Atmospheric ozone and the upper-air conditions," 
Quarterly Journal of the Royal Meteorological Society , v. 79, p. 
39-50, January 1953. 

Ohring, G. and H. S. Muench, "Relationships between ozone and meteoro- 
logical parameters in the lower stratosphere," Journal of Meteor- 
o logy, v . 17, p. 195-206, April 1960. 

Pittock, A. B., "A statistical approach to the synoptic climatology of 
vertical ozone distributions," Journal of Applied Meteorology , v. 

8, p. 308-311, April 1969. 

Prabhakara, C., "Effects of non- photochemica 1 processes on the meridional 
distribution and total amount of ozone in the atrosphere," Month ly 
Weather Review , v. 91, p. 411-431, September 1963. 

Ramanathan, K. R. and R, N. Kulkarni, "Mean meridional distributions 

of ozone in different seasons calculated from umkehr observations 
and probable vertical transport mechanisms," Quarterly Journal of 
the Royal Meteorological Society , v. 86, p. 144-155, April 1960. 

Reed, R. J., "Role of vertical motions in ozone -weather relationships," 
Journal of Meteorology , v. 7, p. 263-267, August 1950. 

Regener, V. H., "On a sensitive method for the recording of atmospheric 
ozone,” Journal of Geophysical Research , v. 65, p. 3975-3977, 
December 1960. 

, "Measurements of atmospheric ozone with the chemiluminescent 

method," Journal of Geophysical Research , v. 69, p. 3795-3800, 

15 September 1964. 

Sticksel, P. R. , The Vertical Distribution of Ozone Over Tallahassee, 

F lor Ida , Ph.D. Thesis, Florida State University, 1966. 



72 



INITIAL DISTRIBUTION LIST 



No. Copies 



1. Defense Documentation Center 20 

Cameron Station 

Alexandria, Virginia 22314 

2. Library 2 

Naval Postgraduate School 

Monterey, California 93940 

3. Department of Meteorology 3 

Naval Postgraduate School 

Monterey, California 93940 

4. Dr. Jerry D. Mahlman 13 



Department of Meteorology 
Naval Postgraduate School 
Monterey, California 93940 

5. Dr. Russell Elsberry 1 

Department of Meteorology 

Naval Postgraduate School 
Monterey, California 93940 

6. Office of the Naval Weather Service 1 

Naval Station (Washington Navy Yard Annex) 

Washington, D. C. 20390 

7. Officer in Charge 1 

Navy Weather Research Facility 

Naval Air Station, Building R-48 
Norfolk, Virginia 23511 

8. Lieutenant Richard H. Stender, USN 1 

C OMNA VMAR I ANAS Box 12 

FPO San Francisco 96630 



9. Professor W. Van der Bijl 1 

Department of Meteorology 
Naval Postgraduate School 
Monterey, California 93940 

10. Professor F. L. Martin 1 

Department of Meteorology 
Naval Postgraduate School 
Monterey, California 93940 



73 



Security Classification 



DOCUMENT CONTROL DATA - R & D 



(Security classification of title, body of abstract and indexing annotation must be entered when the overall report is classified) 



1 originating activity (Corporate author) 


2*. REPORT SECURITY CLASSIFICATION 


Naval Postgraduate School 


Unclassified 


Monterey, California 93940 


26. GROUP 


3 REPORT TITLE 



Numerical Calculations of Stratospheric Ozone Transport 



4 DESCRIPTIVE NOTES (Type of report andjnclusive dates) 

Master’s Thesis; October 1969 

5- AUTHOR(S) (First name, middle initial , last name) 



Richard Henry Stender 



6 REPORT DATE 


7a. TOTAL NO. OF PAGES 


7b. NO. OF REFS 


October 1969 


72 


49 


6a CONTRACT OR GRANT NO. 


9a. ORIGINATOR’S REPORT NUMBER(S) 


6. PROJEC T NO. 






C. 


96. OTHER REPORT NO(S) (Any other numbers that may be aaaigned 
this report) 


d. 







10 DISTRIBUTION STATEMENT 



This document has been approved for public release and sale; its 
distribution is unlimited. 



It. SUPPLEMENTARY NOTES 


12. SPONSORING MILITARY ACTIVITY 




Naval Postgraduate School 




Monterey, California 93940 


13 ABSTRAC T 



The vertical distribution of ozone was observed by the Air Force 
Cambridge Research Laboratories North American ozonesonde network from 
29 April 1963 to 10 May 1963 on a daily basis. This observed data, 
supplemented by twenty-four hour trajectory calculations, was used to 
prepare distributions of ozone mixing ratio on the 200, 100, 50, and 30 
mb surfaces. The ozone distributions were evaluated for the ozone transport 
by mean and eddy motions in the lower stratosphere. The results of these 
calculations show the magnitude of the horizontal and vertical transport of 
ozone into and out of the 30/50, 50/100, and 100/200 mb volume over North 
America. The mean vertical transport, supported by the vertical eddy trans- 
port, is combined to move the ozone primarily downward through the stratos- 
phere. At the tropopause the vertical eddies transport the ozone to the 
troposphere for eventual destruction. 



DD ,'“".“..1473 

S/N 01 01 -807-681 l 



(PAGE 1 ) 



75 



Security Classification 



A-31408 



Security Classification 



1 4 

KEY WORDS 


LINK A 


LINK B 


LINK C 




ROLE 


W T 


ROLE 


W T 


ROLE 


W T 


S tratos phere 
Ozone 

Atmospheric Circulation 
Advect ion 












1 


DD ,"".“.,1473 < "o'* > 










tO|^07-6^| 

7 6 




Security 


Classification 


A 


*3140? 



thesS685 

Numerical calculations of stratospheric 





3 2768 002 02259 2 

DUDLEY KNOX LIBRARY 



